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Abstract 

We present an off-line computational method for 
control allocation design. The control allocation 
function 5 = F(z)r = 6 0 (z) mapping commanded 
body-frame torques to actuator commands is im- 
plicitly specified by trim condition <J 0 (r) and by a 
robust pseudo-inverse problem 

||/ - G(z)F(z ) || < e(z) 

where G(z) is a system Jacobian evaluated at op- 
erating point z, z is an estimate of z, and e(z) < 1 
is a specified error tolerance. The allocation func- 
tion F(z) = ^i(z)Fi is computed using a heuris- 
tic technique for selecting wavelet basis functions 
v and a constrained least-squares criterion for se- 
lecting the allocation matrices F*. The method is 
applied to entry trajectory control allocation for a 
reusable launch vehicle (X-33). 


1. Introduction 

We consider control allocation design for a 
reusable launch vehicle (see Figure 1). The problem 
is described as follows. Let 6 € [R m be a vector of 
actuator conditions (aerosurface positions, engine 
thrust vectoring, etc.) and let z £ IR ri be a vector 
of vehicle operating conditions (angle of attack a, 
mach number M y etc.) such that the the controlled 
vehicle torques 



r may be described as r — / r (<5, z) for some func- 
tion f T . Let r c be a vector of controlled vehicle 
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torques commanded by the attitude control system 
and let z be an estimate of vehicle operating condi- 
tions z. The control allocation problem is to design 
a control allocation law 5 c (z,r c ) specifying actua- 
tor commands 5 C such that vehicle body torques 
f T {z y (5 c (z, t c )) match the commanded torques r c to 
within a specified tolerance. 

We present an off-line computational approach 
for a gain scheduled linear affine control allocation 
law 5 c (z, r c ) = F(z)r c + 6o{z) where do(*) is a trim 
function selected such that / r (do(z),z) = 0 and 

Tl w 

F(z) = Y,Ms)Fi (1.1) 

1=1 

is a gain scheduled control allocation matrix com- 
puted from basis (wavelet) functions 'F = 
and their corresponding weighting matrices T — 

(W=Y 


2. Control Allocation Design 

Let H be a region of operating points (z £ fi) 
and define 

G{z)£VsM6,z)\ i=i0(s) :ze n 

Ideally, the control allocation matrix F(z) is se- 
lected to be a pseudo-inverse of G(z), 

G(z)F(z) = /. 

However, lack of precise measurements of operat- 
ing conditions z and actuator saturation issues can 
prevent achievement of this ideal condition. 

The use of an estimate z for gain scheduling re- 
quires that we relax the above pseudo-inverse condi- 
tion. Let e r (z) be a specified error tolerance defined 
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on z £ U and let Q — Q r > 0 be a selected matrix 
to define a quadratic norm \\z\\ Q = 1 y/z T Qz, Then 
a revised control allocation performance condition 
may be written as 

\\l-G(z)F(z)\\<e r (z):\\z-z\\ Q <e z ( 2 . 1 ) 

We require tr(z) < 1 since this condition guaran- 
tees t t t c > 0, i.e., the resulting controlled vehicle 
torques have “the correct sign.” which is a neces- 
sary condition for several attitude control strategies 

([5], [6]). 

Actuator saturation issues may be dealt with in 
part by incorporating the constraint 

\\W s F{z)W r \\ < l,z £ ft (2.2) 

where W 6 and W r are weighting matrices selected 
based on actuator saturation limits and anticipated 
maximum torque command magnitudes, respec- 
tively. Notice that equation (2.2) is a convex con- 
straint on the linear-affine matrix function F(z). 

In the absence of constraints (2.1) and (2.2), the 
ideal control allocation problem may be relaxed to 
a simple least-squares optimization 

min [I \\I - G(z)F(z + h)\f dhdz 

Jzen y/i:||/i|| g <€, 

^ ( 2 - 3 ) 
The use of constraints (2.1) and (2.2) yields a con- 
vex programming problem over z,z. 

We propose a candidate (off-line) solution pro- 
cedure that iteratively selects a set of basis wavelet 
functions # and a corresponding set of weighting 
matrices T. Our candidate procedure assumes that 
the optimization (2.3) subject to constraints (2.1) 
and (2.2) is feasible; that is, there exist V L, T that 
satisfy the constraint conditions. The procedure is 
as follows. 

Algorithm 2.1 Control allocation design proce- 
dure. 

Inputs Target error tolerance function e r (r), sys- 
tem Jacobians (7(.c), measurement uncer- 
tainty parameters Q = Q r > 0 and e, > 0. 

Outputs Wavelet set ^ and a corresponding set of 
weighting coefficient matrices T to specify the 
allocation function F(z) per equation (1.1). 


1 . *:= -1 
2. do 

(a) £:=£+! 

(b) if A- = 0 then select an initial set of can- 
didate wavelets 'I'o* 

(c) else Update to obtain a new set of 

candidate wavelets 

(d) end if 

(e) Compute a set of matrices T k corre- 
sponding to satisfying the minimiza- 
tion (2.3). 

(f) Define F k (z) from equation (1.1) for ^ £ 

F t £ T k 

3. until F k satisfies constraints (2.1) and (2.2). 
Remark 2.1 Let 

e* = max||7 — G{z)F(z)\\ - e r (z) 

: z,z £Sl,\\z-z\\ Q <€ z 

The wavelet sets should be selected such that 
the sequence e k is monotone decreasing. In princi- 
ple, one may select the set of wavelets $ k so that 
they provide increased resolution at values of z £ H 
where the constraint (2.1) is violated by F k -i(z) 
corresponding to $*_i, By consequence, the 

updated wavelets provide additional degrees of free- 
dom to the unconstrained optimization (2.3) so that 
the computed values T k yield a function F(z) that 
is “nearer” to a feasible solution to the constrained 
optimization. 


3. Computational approach 

We now consider Algorithm 2.1 in further detail. 
Since the function f T (z y S) is often available only in 
sampled form through, e.g., wind tunnel data, w'e 
select a set of N operating points Z = £ f2 

in the operating region Q. Associate with Z the 
set of Jacobians Q = {G, = G(z l ) y z t £ Z}. W ; ith 
the function G(z) thus sampled, given a set of n u , 
wavelets ^ the unconstrained control 

allocation problem (2.3) becomes a matrix-valued 
(dense) least squares problem of the form 

mm\\AX-B\\ 
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whore 


F x ' 


■ / ' 


and B = 




i _ 


and each block-row .4,. of .4 is of the form 

[ ••• rp n „(zj)Gi] (3.1) 

where \\zj - Zi\\ Q < e z , which may be solved read- 
ily through a QF-factorization [1]. Further, if the 
dimension of the minimization is not too large, re- 
cently developed convex programming techniques 
[3], [9] allow the norm constraint 

||/-C(z i )F(£ i )|| < € T {Z) 

A'tjZ* £ 2, 1 1 z j Zj || ^ 6 - 

to be directly included in Algorithm 2.1 step 2e. 
Should the dimension of the minimization (3.1) be 
too large for practical direct application of convex 
programming techniques, one may employ an iter- 
ative algorithm to update individually the matrix 
values Fi £ J r k) holding the other unknowns fixed. 
More precisely, we may replace step 2e of Algorithm 
2.1 with the following procedure. 

Algorithm 3.1 Tune weighting matrices 

1. Compute T k from the unconstrained least 
squares minimization (3.1). 

2. for each wavelet i/jj k ^ £ 

(a) Identify sample points z 3 £ H contained 
in the support of xi)\ k \ 

(b) Select weights ci,C 2 >0 and update the 
weighting matrix F [ x k) £ T k to solve 

min + 

subject to || / - G(zj)F(zj) || 

< t r (z) + g 
\\W 6 F{z 3 )W t \\ < 1 + t 

3. end for 

Since the wavelets £ V are assumed to be 
of finite support, for wavelets of higher degrees of 


resolution this reduces not only the column dimen- 
sion of A but also its row’ dimension, and allow’s 
the incorporation of the norm constraint (see equa- 
tion (2.1)) without undue increase of computational 
complexity. 

Remark 3.1 Notice that Step 2b above does not 
directly enforce the constraints (2.1) and (2.2), but 
rather induces a trade-off betw r een them in the com- 
putation. This is done to ensure feasibility of the 
posed convex optimization problem. Algorithm be- 
havior may be tuned as needed by adjustment of 
selected coefficients cj , co or by adding other con- 
vex constraints, e.g., g > 0, as needed. 

3.1. Wavelet decomposition 

Formal development of a wavelet decomposition 
of a general function F(z) [7] involves w’avelet nor- 
malization and the definition of a multidimensional 
filter bank. Because the function F(z) is only im- 
plicitly defined by the error criteria (2.1) and (2.3), 
the structure of the control allocation problem does 
not readily allow exploitation of wavelet properties 
of (bi-)orthogonality, etc., nor does does it admit di- 
rect application of a subset selection algorithm [2]. 
Reeves’ Backward Greedy Algorithm [4] may be ap- 
plied to the unconstrained problem (3.1), but this 
approach does not make use of the finite-support 
property of the set of wavelets nor does it read- 
ily admit inclusion of the norm constraint (2.1) re- 
stricted to the set of operating points 2 . 

For our purposes, it suffices to develop a heuris- 
tic algorithm to identify a small set of wavelet ba- 
sis functions {x such that the error tolerance 
(2.1) can be achieved. For clarity in exposition, we 
shall use tw’o-dimensional wavelets y) w'here 
x 

= Tz is a vector drawm from a 2-D hvper- 

plane in IR n . (Extension of our technique to higher 
dimensions is not difficult, but rapidly increases the 
computational complexity of the problem in terms 
of both off-line computation time and required stor- 
age in flight software.) We select the set of two- 
dimensional wavelets 

xp {%jkl) (x, y) = ruin (x), v [kt) (g) j 

where i,j, k , f are integers, 

^ (,j) (x) = V{2~i(z - j)). 
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and v(-) is the mother wavelet function (see Figure 

2 ) 

( x+l x€ [-1,0] 

ti>(x) : IR ->• IR : x -» ^ 1-x x 6 [0, 1] 

( 0 else 

(Notice that we do not normalize ||^(*J*0|| = l ; 
justification for this will be seen Algorithm 3.2 in 
the next subsection.) This class of wavelet ba- 
sis functions may alternatively be interpreted as a 
Sugeno-Takagi fuzzy logic (STFL) function or as a 
first-order finite element solution; the use of convex 
programming in STFL based dynamic control for a 
fixed set of basis functions 4> is addressed in [8] . We 
will occasionally write F(x,y) rather than F(z). 

3.2. Selection of wavelets 

Since the matrix values T = are easily 

computed as the solution of a linear least-squares 
problem, we consider now the selection of the set of 
basis wavelet functions and its cardinality 

n w . We select the sequence of wavelet sets 'F* in 
Algorithm 2.1 in an iterative process of subdividing 
two-dimensional wavelets y) with low in- 

dex values i and k into ‘‘higher frequency” wavelets. 
We shall make use of the following terms: 

Definition 3.1 Consider an arbitrary wavelet 
^(x,y) = xp (ijke) {x,y). 

1. The center of ip is the point 

(£,$) = (j 2- { ,f2-*) 

at which the wavelet has its maximum value. 

2. The x axis children of t/>(x,y) are the wave- 
lets •■>'•*’*>} with j' = 2 j - 1 , 2j, 2j + 1. 

Notice that the support of the wavelet ip is the 
union of the support of its x-axis children. 
Similarly, the j / — axis children of w(x, y) are 
{^(*o,*Ti,r')} f = 2f- l,2f,2f + 1. The 
support of the wavelet ip is the union of the 
support of its y- axis children. 

3. We say that a wavelet ip [ 'i kt) {x,y) has uni- 
form resolution in both axes if i = k. If v has 
uniform resolution we say that its resolution 
level is i. 

4. Suppose that y> h;is uniform resolution. Its 
uniform resolution child wavelets, or simply 


its children, are the x-axis children of its 
y-axis children (or, equivalently, the y-axis 
children of its x — eixis children) whose centers 
are contained in the region ft. Each wavelet 
has up to nine children. 

5. Consider a set of wavelets 4» with uniform res- 
olution. The wavelet ip with resolution level 
i is redundant with 4* if (1) 0 g <F, or (2) its 
children are all redundant with 4'. Alterna- 
tively, ip is redundant with 4' if the support of 
ip can be written as the union of the supports 
of wavelets in a set 4 r ' C <F whose resolution 
level exceeds i. 

We now present a refined version of Algorithm 2.1. 

Consider the N sample points Z = {z J }N =1 g ft. 

Assume without loss of generality that ft = [0 11 x 

[o,i]. 

Algorithm 3.2 Control allocation design proce- 
dure 

1. k = 0; tfo is selected to be the 2 P = 4 basis 
wavelets with resolution level 0 centered at 
the vertices of ft, 

= (W 0 ) (x,y). ...,t^ 0 ) (x,y)} 

_ |^( 0 , 0 , 0 , 0 ) ^( 0 , 0 , 0 , 1 ) , ( 0 , 1 , 0 , 0 ) 

^( 0 , 1 , 0 , 1 )}. 

2. Compute a set of matrices T§ corresponding 
to »Fo satisfying the minimization (3.1). 

3. Define F 0 (z) from equation (1.1) for t />,• g 4 > 0 
Fi 6 F 0 

4. while Ft does not satisfy the constraint 

||/-G(z i )F(£ l )|| < e T (z) 



(a) Compute the error surface 

e*fy.) = max||/ - Gfy I )F(£)|| 

— £t(z ,) 

: Zi,z, g Z, ||r, - ;,|| < <. 

and define e k (z,) = max{0, e*(z,)}. 
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(b) Identify a wavelet v E V k with maximal 
inner product with the error surface e k . 

(c) k := k + l 

(d) Construct V k from by replacing xp 
with its children. 

(e) Define F k (c) from equation (1.1) for \p x E 

Fi E T k 

5. end while 


Remark 3.2 Notice that e* is nonnegative at all 
Zi £ Z and is positive only at points z t € Z where 
the norm constraint (2.1) is not met. {^ ijkt \e k ) = 
0 => the constraint (2.1) is met at all points z { in 
the support of xp^ k( ^ 


Remark 3.3 Since the wavelets are not 

normalized to have unit energy, the inner product 
criterion in step 4b is biased toward subdividing 
low-resolution wavelets. This heuristic process of 
subdivision quickly increases the level of resolution 
in those areas of Cl where it is needed while using 
low resolution levels where possible, reducing the 
computational requirements of the final control al- 
location flight code. 


4. Design Example 

We illustrate our design procedure in applica- 
tion to the X-33 experimental technology vehicle 
entry trajectory control allocation design. From 
the X-33 main engine cut off (MECO) to its tar- 
get area energy management (TAEM) interface, the 
X-33 actuators include two pairs of elevons (left, 
right; outboard, inboard) and two body flaps (left, 
right). The left and right elevon pairs operate as 
a unit so that there are four deflection commands 
(left elevon, right elevon, left body flap, right body 
flap), 6 C E IR 4 . Due to details of the vehicle, the 
trim deflection function for the X - 33 is selected 
as 


^left elevon 


" <5, • 

^right elevon 


tie 

^left flap 


M:) 

bright flap 


Ms) 


where 8 e is a (constant) elevon bias position and 
ft(z) is a scalar- valued flap pitch trim function. 


(Roll and yaw torques resulting from <S 0 (z) are all 
zero.) A set of eight reaction control system (RCS) 
thrusters are also available for entry maneuvers; for 
simplicity (and due to details of RCS thruster op- 
eration) these are neglected in this study. 

The control allocation is scheduled on the an- 
gle of attack a and the Mach number M. 3-DOF 
simulations and dispersion (Monte Carlo) runs indi- 
cate that the range of these values is a E [10°,45°] 
and M E [3, 10]. For the purposes of wavelet de- 
composition these ranges were normalized to fl = 
[0, 1] x [0, 1]. A set of 1024 operating points Zj E Cl 
were selected by superimposing an evenly spaced 
32 x 32 rectangular mesh on the above intervals. A 
plot of condition k(G(z)) of the Jacobian G(z) from 
the surface deflections 6 to body torques r at trim 
conditions is shown in Figure 3. Based on the data 
in Figure 3 and other analysis of achievable per- 
formance, the target tolerance e r (z) was selected 
as shown in Figure as shown in Figure 4 with val- 
ues between 0.5 and 0.71. The initial thresholded 
error surface e\ for n w = 4 is shown in Figure Fig- 
ure 5. Two iterations later ( k = 3) there are a 
total of n w = 9 wavelets in ^ 3 and the error sur- 
face e 3 is as shown in Figure 6; notice the decrease 
in maximum error above threshold and in the por- 
tion of Cl that meets the error criterion (2.1) (points 
at which e 3 =0). The target tolerance e r (z) wa s 
achieved after a total of four iterations with a to- 
tal of n w = 14 wavelets in the database; a plot 
of ma Xjipj 4 \z) for xv^ E ^4 is shown in Figure 
7; peaks in the plot correspond* to the centers of 
wavelet functions. (Use of least-squares optimiza- 
tion without re-tuning individual weights as in Al- 
gorithm 3.1 required n u , = 53 wavelets to achieve 
similar performance.) 

The computed control allocation function F(z) 
was tested in open-loop simulation by command- 
ing a sawtooth wave in roll, pitch, and yaw torques 
(see Figures 8- 10). Torque command channels are 
coupled, but have acceptable performance as shown 
in the Figures. Aerosurface deflections (with ideal 
pitch trim function Jo(^)) are shown in Figure 11. 


5. Conclusions 

The design of control allocation for tail-less, lift- 
ing body reusable launch vehicles requires both ro- 
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bust ness to measurement and modeling uncertainty 
and a computationally feasible implementation. We 
propose the off-line computation of a wavelet-based 
decomposition of an implicitly specified control al- 
location function 

F( x ) — rJji(z)F l + (5 0 (c) 

i 

We select piecewise linear wavelets through an it- 
erative subdivision process in order to meet a max- 
imum error criterion. The matrix weights {FJ 
are computed through the solution of a linear least 
squares problem. The method was applied to entry 
trajectory control allocation for the X-33 technol- 
ogy testbed and appears promising. 
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Figure 1: Block diagram of control allocation prob- 
lem. It is desired to develop a set of con- 
trol deflection commands S c such that the 
resulting body torques r are “near” to the 
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Figure 2: 1-D piecewise linear mother wavelet. 
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Figure 3: Pointwise condition number of the de- 
flection to body torque Jacobian matrix 
G(a, M). Omitted values rapidly rise to a 
condition number k ( G(a , A/)) a 160. Cor- 
responding values of G show less than 1% 
of yaw authority than that available in the 
pitch and roll channels. Mean and standard 
deviation of the condition number values 
shown in the plot are both approximately 
15. 


Figure 5: Initial error thresholded by target tolerance 
value 
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Figure 4: Maximum error tolerance for control allo- 
cation matrix computation 


Figure 6: Threshold error surface after two iterations 
of wavelet subdivision 
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Wavelet basis [09-Mar-2000 13:15:40] 



Figure 7: ma xjipj(z) for wavelet basis functions selected for control allocation function F( 



Figure 8: Control allocation test in simulation: roll command torque r c and body torque r. 
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Figure 9: pitchtorque 












